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00 , In this paper, we study two important metrics in multiple-input multiple-output (MIMO) time-varying Rayleigh 

flat fading channels. One is the eigen-mode, and the other is the instantaneous mutual information (IMI). Their second- 
order statistics, such as the correlation coefficient, level crossing rate (LCR), and average fade/outage duration, are 
jyT ' investigated, assuming a general nonisotropic scattering environment. Exact closed-form expressions are derived and 

^, ' Monte Carlo simulations are provided to verify the accuracy of the analytical results. For the eigen-modes, we found 

they tend to be spatio-temporally uncorrelated in large MIMO systems. For the IMI, the results show that its correlation 
►^ ■ coefficient can be well approximated by the squared amplitude of the correlation coefficient of the channel, under 

CO , certain conditions. Moreover, we also found the LCR of IMI is much more sensitive to the scattering environment 

^— s , than that of each eigen-mode. 
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I. Introduction 

The utilization of antenna arrays at the base station (BS) and the mobile station (MS) in a wireless communication 
system increases the capacity linearly with min(A^, iVp), under certain conditions, where Nt and A'r are numbers 
of transmit and receive antenna elements, respectively, provided that the environment is sufficiently rich in multi-path 
components [1][2]. This is due to the fact that a multiple-input multiple-output (MIMO) channel can be decomposed 
to several parallel single-input single-output (SISO) channels, called e/gen-channels or e/ge«-modes\ via singular 
value decomposition (SVD) [2]-[9]. 
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For a SISO channel, or any subchannel^ of a MIMO system, there are numerous studies on key second-order 
statistics such as correlation, level crossing rate (LCR), and average fade duration (AFD) [10]-[14]. However, to 
the best of our knowledge, no such study on the eigen-channels of a MIMO system is reported in the literature, 
possibly due to the lack of knowledge regarding the joint probability density function (PDF) of e/gen-channels. 

Regarding another important quantity, the instantaneous mutual information (IMI), only some first-order statistics 
such as the mean, variance, outage probability and PDF are studied [9][15]-[18]. Clearly, those statistics do not show 
the dynamic temporal behavior, such as correlations, LCR and average outage durations (AOD) of the IMI in time- 
varying fading channels. It is known that IMI can be feedbacked to the rate scheduler in multi-user communication 
environments, to increase the system throughput [16], where only the perfect feedback is considered. However, it is 
hard to obtain perfect feedback in practice due to the time-varying nature of the channel, which makes the feedbacked 
IMI outdated. In this case, the temporal correlation of IMI can be used to analyze the scheduling performance with 
outdated IMI feedbacks. Furthermore, one can improve the rate scheduling algorithm by exploring the temporal 
correlation of IMI. 

Several second-order statistics such as the correlation coefficient, LCR and AOD of IMI in single-input single- 
output (SISO) systems are reported in [19] and [20]. For MRC-like MIMO systems, they are investigated in [20]. 
However, there are a limited number of results for a general MIMO channel. In [21], some simulation results 
regarding the correlation coefficient, LCR and AOD are reported, without analytical derivations. In [22], lower and 
upper bounds, as well as some approximations for the correlation coefficient of IMI are derived, without exact 
results at high SNR. A large gap between the lower and upper bounds and large approximation errors are observed 
in [22, Figs. 2, 5]. 

In this paper, we extend the results of [20] to the general MIMO case, using the joint PDF of the eigenvalues 
[23]. Specifically, a number of second-order statistics such as the autocorrelation function (ACF), the correlation 
coefficient, LCR and AFD/AOD^ of the e/ge«-channels and the IMI are studied in MIMO time-varying Rayleigh 
flat fading channels. We assume all the subchannels are spatially independent and identically distributed (i.i.d.), 
with the same temporal correlation coefficient, considering general nonisotropic scattering propagation environments. 
Closed-form expressions are derived, and Monte Carlo simulations are provided to verify the accuracy of our closed- 
form expressions. The simulation and analytical results show that the eigen-modes tend to be spatio-temporally 
uncorrected in large MIMO systems, and the correlation coefficient of the IMI can be well approximated by the 
squared amplitude of the correlation coefficient of the channel if |iVr — Nii\ is much larger than miii (iVr, Nn). In 
addition, we also observed that the LCR of IMI is much more sensitive to the scattering environment than that of 
each eigen-mode. 

The rest of this paper is organized as follows. Section II introduces the channel model, as well as the angle- 
of-arrival (AoA) model, figen-channels of a MIMO system are discussed in Section III, where Subsection III-A 



^In this paper, each subchannel represents the radio link between each transmit/receive pair of antennas. 
'Note that AFD is used for eigen-channels. whereas AOD is used for MIMO IMI. 



is devoted to the derivation of the normaUzed ACF (NACF) and the correlation coefficient of e/ge«-channels of a 
MIMO system, whereas Subsection III-B focuses on the LCR and AFD of the e/ge«-channels. The MIMO IMI 
is investigated in Section IV, in which Subsection IV-A addresses the NACF and the correlation coefficient of 
the MIMO IMI as well as their low- and high-SNR approximations, whereas Subsection IV-B studies the LCR 
and AOD of the MIMO IMI using the well-known Gaussian approximation. Numerical results and discussions are 
presented in Section V, and concluding remarks are given in Section VI. 

Notation: -^ is reserved for matrix Hermitian, * for complex conjugate, j for V^-T, E[-] for mathematical 
expectation, !„ for the m x m identity matrix, || • \\p for the Frobenius norm, 5ft[-] and 3[-] for the real and 
imaginary parts of a complex number, respectively, and /^(x) for [f{x)] . Finally, ie[m, n] implies that t, m and 
n are integers such that m < t < n with m < n. 

II. Channel Model 

In this paper, an Nji x Nt MIMO time-varying Rayleigh flat fading channel is considered. Similar to [15], we 
consider a piecewise constant approximation for the continuous-time MIMO fading channel matrix coefficient H(i), 
represented by {ll{lTs)} [^^, where Tg is the symbol duration and L is the number of samples. In the sequel, we 
drop Ts to simplify the notation. In the l^^ symbol duration, the matrix of the channel coefficients is given by 



H(0 



hNR,l{l) ■■■ ^Affl.AfrlO 



Je[i,L]. (1) 



We assume all the NtNh subchannels {/i„^ „j(/),/ e [lii]}(„"Li n =i) ^^^ i.i.d., with the same temporal corre- 
lation coefficient, i.e., 

E[hmn{l)hpg{l - i)] = Sm,pS„^qPhii), (2) 

where the Kronecker delta 5m,p is 1 or when m = p or m ^ p, respectively, and Phii) is defined and derived at 
the end of this section, eq. (4). 

In flat Rayleigh fading channels, each hn^.n^il), I G [1, L], is a zero-mean complex Gaussian random process. In 
the ^*'' interval, /i„^,nj(^) can be represented as [13] 

(3) 

= anr.,nt (0 exp[-j$„^,„j (0], 
where the zero-mean real Gaussian random processes /i^j ,14(0 ^"^ ^n^ nti^) ^^ f^he real and imaginary parts of 
^nr,nf(0' respectively. a„^.„j(Z) is the envelope of /i„,..„j(/) and ^nr.nt{l) is the phase of ft,„^,„j(/). For each /, 
ttrir, nt(0 has a Rayleigh distribution and $„^,„j(Z) is distributed uniformly over [— tt, tt). Without loss of generality, 
we assume each subchannel has unit power, i.e., IE[a^ „ (l)] = 1. 

Using empiricaUy-verified [13] multiple von Mises PDF's [19, (4)] for the AoA at the receiver in nonisotropic 
scattering environments, shown as Fig. 1 of [19], the channel correlation coefficient of hn^^nA^)^ ^nrjUt, is given 



by [19, (7)] 



Pn{l)^nhn.,nAl)K^,n,{l-^)]. 

N 



Io(V< ~ 47r2/^i22;2 + jAn^,„fmTs cos 0„) (4) 

where Ik{z) = 7/9 e^'^"'"" cos{kw)d'w is the fc* order modified Bessel function of the first kind, 0„ is the mean 
AoA of the n* cluster of scatterers, k„ controls the width of the n* cluster of scatterers, P„ represents the 
contribution of the n* cluster of scatterers such that X]n=i ^n = li < P„ < 1, A' is the number of clusters of 
scatterers, and fn is the maximum Doppler frequency. When k„ = 0, Vn, which corresponds to isotropic scattering, 
(4) reduces to phii) = loiftT^ loiTs) = Jo{2Tr foiTg), which is the Clarke's correlation model. 

III. Eigen-Channels in MIMO Systems 

We set AI = inm{NT,Nii) and N = max(i\^, iVp). Based on singular value decomposition (SVD) [2]-[9], 
H(^) in (1) can be diagonalized in the following form 

H(0 = U(OS(OVt(0, (5) 

where Y{1), whose dimension is Nt x M, satisfies V^(/)V(/) == Im, U(/), which is NrxM, satisfies U'f(/)U(/) = 
Ia/, and S(/) is a diagonal matrix, given by S(/) — diag [si (/),•• • ,sm(0]' ii^ which Sm{l), m G [Ij-^j'^] is the 
TO*^ non-zero singular value of H(/). 

We define Xm{l) = s^(/), Vm. Therefore Xm{l) is the m'^ non-zero eigenvalue of H(^)Ht(/). We further 
consider {Am(0}m=i ^^ unordered non-zero eigenvalues of H(/)Ht(/). Therefore, the MIMO channel H(/) is 
decomposed to M identically distributed e/ge«-channels, {A„i(^),^ G [Ij-^llm^i' by SVD, as shown in Fig. 1. For 
M = 1, there is only one e/§e« -channel, which corresponds to the maximal ratio transmitter (MRT) if Nji = 1, or 
the maximal ratio combiner (MRC) if Nt — 1. In each case, we have N i.i.d complex Gaussian branches. 

Since all the e/ge«-channels have identical statistics, we only study one of them and denote it as X{1), I E [1, L]. 
To simplify the notation, we use X and Y to denote A(/) and X{1 — i), respectively. The joint PDF of X and Y is 
given in (6) of [23], 






kill 



M\l~gi)e'^ ^^ {k + u)lQf M2 ^A^^^\{^k + u)l{l + v)l 

X [[Ll{x)L'{{y)f + [mx)Ll{y)f - [gf-""^ + ^f "'^J Ll{x)L'{{x)Ll{y)L'({y)]Y (6) 

where Lf^{x) = ^e^x^"^^ (e^^x"+") is the associated Laguerre polynomial of order n [24, pp. 1061, 8.970.1], 
u = N — M, and Qi = \ph{i)\, where ph[i) is given in (4). The joint PDF in (6) is very general and includes many 
existing PDF's as special cases [23]. 

• By integration over y, (6) reduces to the marginal PDF 

A/-1 , 

m=o ^ ' 



which is the same as the PDF presented in [2]. When AI — 1, (7) further reduces to 

P{-) = (^-"-^e-, (8) 

which is the x^ distribution with 2N degrees of freedom [25, (2.32)], used for characterizing the PDF of 
outputs of MRT or MRC [26]. 
With M = 1, (6) reduces to 



(-»)^exp(-|ij)/„_,(Jf^) 



which is the joint PDF of outputs of MRT or MRC at the ^*'' and {I - i)*'^ symbol durations [27]. It includes 
(3.14) of [25] as a special case'*. Furthermore, when A^ = 1, i.e., a SISO channel, (9) simplifies to 

1 / x + y\ (2g^^/xy" 

2 '^^P ^"i 2 ^0 -; 2 

^-gf \ i-gfj \ ^~gf 

which is identical to (8-103) [28, pp. 163], after a one-to-one nonlinear mapping. 
In the following subsections, we study the normalized correlation and correlation coefficient of any two eigen- 
channels, defined by, respectively. 



p{x,y)^-^ — 12 '^^Pl^-; — 3 Mo ( / „2 ) ' (10) 



and 



''"■"^'^ " yiElAlMv/EMl-)!' ^ ^ 

E[A™(OA„G-z)]-E[A,„(/)]E[A„(/-z)] 

yOm.ri(«) = , , , (12) 

y'E [A2„(0] - {E [\m{l)]f^J^ [\l{l - i)] - {E [A„(/ - l)]f 

A. Normalized Correlation and Correlation Coefficient of Eigen-Channels 

To derive the normalized correlation and correlation coefficient between any two eigen-channels, we need the 
following lemmas. 

Lemma 1: The first and second moments of the tti"^ eigen-channel are respectively given by 

E[A„,(0]=iV, (13) 

E[A^(0]=iV(iV + M). (14) 

Proof: See Appendix I. D 

Lemma 2: The autocorrelation of the m"' eigen-channel, defined as r^n.-mii) — E [Am(0'^m(^ — «)], is given by 

rra.mil) =N^ + ^, * ^ 0. (15) 

Proof: See Appendix II D 
Lemma 3: The cross-correlation between the m*'' and n"' eigen-channels, defined as r™ „(i) = E [Am(/)A„(/ — i)]. 



is given by 



N^ - N, i = 0, 
rni.n{i) = { , n^ m. (16) 



Eq. (3.14) in [25] is developed for real uncorrelated Gaussian random variables. 



Proof: See Appendix III. D 

Based on Lemmas 1-3, we obtain the closed-form expressions for (1 1) and (12), which are given in the following 

theorem. 

Theorem 1: The normalized cross-correlation and the correlation coefficient between wy^ and n*'' eigen-channels, 

defined in (11) and (12), are respectively given by 

J\/~(M+l)(l~^„.„) „ 

N+M ' * — "7 



MN+gj 
MN+M^ ' 



i^O, 



and 



Af+l 



PmA^)={ , '' ^' ''""^' '~"' (18) 

Proof: From Lemma 1, it is straightforward to see that the eigen-channel is stationary in the wide sense. More- 
over, all the eigen-channels have the same statistics, therefore we have E [A„(Z — i)] = E [Am(0] ^"d E [A^j(/ — i)] = 
E [A^„(Z)] , Vm, n e [1, M] and Wl, i. By plugging (14)-(15) into (11), we obtain (17). Finally, substitution of (13)- 
(15) into (12) results in (18). D 

From (17) and (18), we have the following interesting observations. 
• If M is greater than 1, the normalized correlation and the correlation coefficient are not continuous at z = 0, 

as f,„,„(l) and Pm,„(l) do not converge to rm,„(0) = p,n.n{^) = 1 as 2; -^ 0, Vm, n. 
m If M is large, all the M eigen-channels tend to be spatio-temporally uncorrected, due to 

lim Pm.n{i) = Sjn,nSi.Q- (19) 

M— *oo 

As an example, with isotropic scattering, (17) and (18), respectively, reduce to 

M-(M+l){l-S^_^) „ 

rm,nW = S , (20) 

MN+jS{2^fm%) . , „ 

MN+AP ' '7^"' 



and 



1 - ^ (1 - <5,„,„) , * = o, 



Pm,n{i) = S , „ (21) 

'■2ttJi/ 



J;^{27rfotTs) .^Q^ 



B. LCR and AFD of an Eigen-Channel 

In this subsection, we calculate the LCR and AFD of an eigen-channel at a given level. To simplify the notation, 
the eigen-channel index rn is dropped in this subsection, as the derived LCR and AFD results hold for any eigen- 
channel. 

1) LCR of an Eigen-Channel: Similar to the calculation of zero crossing rate in discrete time [29, Ch. 4], we 
define the binary sequence {Zi}i^^, based on the e/gen-channel samples {X{l)}i^^, as 

|l, ifA(/)>A,h, 
Zi = I (22) 

0, ifA(/)<Ath, 



where Ath is a fixed threshold. The number of crossings of {)^{l)}i^i with Ath, within the time interval Tg < t < LTg, 
denoted by Dau,, can be defined in terms of {Zi}fL^ [29, (4.1)] 



1=2 

which includes both up- and down-crossings. 

After some simple manipulations, the expected crossing rate at the level Ath can be written as 

E[Z?aJ 2Pr{Zi = 1} - 2Pr{Zi = 1, Z;_i = 1} 



(23) 



(24) 



where P,!} is the probability of an event. Therefore, the expected down crossing rate at Ath, denoted by N\{Xih 
is half of (24), given by 



NxiX±) 



0A(Ath) - ^\i\h) 



(25) 



where 4'\{Xth) = Pr{Zi = 1} and iy9A(Ath) = Pr{Zi = l,^;-i = 1}. Analytical expressions for (/)A(Ath) and 
f\{^t\\) are stated in the following theorem. 

Theorem 2: For a given threshold Ath, 0a (Ath) and 'p\{\t\\) are, respectively, given by 






"-' " " ^!(:::^pcrpr(p+g+.+i,Ath) 



m— Op— Og— 



and 



oo A/-1 



'fixi^h) = (f>l{>^th) + 



\,n 2^ 2^ (n 



Jlklgl'^'-"^ 



j=M fc=0 ^■' ' ^ ' 



EE :f:r^.tl! r(p+g+.+i,Ath) 

,p=0 q-- 



(26) 



(27) 



where r(a, z) = J, i"^^e^*dt [24, pp. 949, 8.350.2] is the upper incomplete gamma function, (") is the binomial 
coefficient, given by ^i-^^^s, , and q\ = |p/i(l)|, defined before, i.e.. 



Q\ = 



E^»^ 

n=l 



/o('«n) 

Proof: L'^{x) is a polynomial of order n, and can be represented as [24, pp. 1061, 8.970.1] 



K{^)^Y. 



k=0 



-k 



k\ 



(28) 



(29) 



By plugging (29) into (7), the univariate PDF of an eigen-channel, and integrating over x from Ath to oo, we obtain 
(26). Similarly, substitution of (29) into (64), the bivariate PDF of an eigen-channel, and integration over x from 
Ath to oo results in (27). D 

By plugging (26) and (27) into (25), we obtain the expected crossing rate at the level Ath. 

2) AFD of an Ei gen- Channel: The cumulative distribution function (CDF) of X{l),yi, is obtained as 



PA(Ath) ^Pr{X< Ath} = 1 - 0A(Ath), 

where 0a (Ath) is given in (26). 

The AFD of the e/ge«-channel {X{l)}i^i is therefore given by 

PA(Ath) _ [l-0A(Ath)]T, 



(30) 



tA(Ath) = 



^A(Ath) 0A(Ath)-'/'A(Ath)' 



(31) 



where 4>\{^th) and (p\{Xth) are given in (26) and (27), respectively. 



IV. MIMO IMI 

In this section, the NACF, the correlation coefficient, LCR and AOD of IMI in a MIMO system are investigated 
in detail. In the presence of the additive white Gaussian noise, if perfect channel state information {H(/)}j^j^, is 
available at the receiver only, the ergodic channel capacity is given by [2] [9] 



C = E 



Indet ( Iat^ + ^HiH/ 



Nt 



(32) 



in nats/s/Hz, where 77 is the average SNR at each receive antenna, and Hj denotes H(/). 

In the above equation, at any given time index I, In dot [iNn + n~^I-'^i) ^^ ^ random variable as it depends 
on the random channel matrix H/. Therefore 

V 



Indet I 



Wb 



Nt 



-HzHt 



^ = 1,2, 



(33) 



is a discrete-time random process with the ergodic capacity as its mean. 

By plugging (5) into (33), we can express the IMI in terms of A/ eigenvalues as 



M , 

^hih+^A™(0), z = i,2,. 

m — 1 ^ 



(34) 



A. NACF and Correlation Coefficient of MIMO IMI 

In this subsection, we derive exact closed-form expression fors the NACF and the correlation coefficient of MIMO 
IMI, and their approximations at low- and high-SNR regimes, using the following lemmas. 
Lemma 4: The mean and second moment of Ii are respectively given by (35) and (36) 



M-l 



m— p— g— 



\m—p/ \m — q/ 



{■m+iyy.p\ql{-l) 



p+q ^,J\ jj 



0,1 

0,0,p + q + i' + l. 



(35) 



A/-1 






ni=0 p=0 q=0 j=0 ( jj- 



p+q+v+l 

— I {m + i')\plql 



'3,4 



M-1A/-1 

EE 



j\k\ 



(j + j.)!(fc + z.)! 



k-q) 



'-;){ 



^ ^ plql(-l)P+i 



G. 



;,i Nt 



j=0 k=0 

where G is Meijer's G function [24, pp. 1096, 9.301]. 
Proof: See Appendix IV. 
Lemma 5: The ACF of MIMO IMI, defined as rz{i) = E[X;X;_i], is shown to be 

J k /j+i'-) /k+u 




V 

0,1 

0,0,p+q+iy + l 



00 M-l ... , 2(j-k) 



j=M k=a w I / ^ i ; 
Proof: By plugging (29) into (64), and using (72), we obtain (37) immediately 



EE 



k-ql ^3,1/ ^T 



p^o^t^pi'^K-ir^ '-'{v 



0,1 

0,0,p+<7+i^+ly 



(36) 

D 

(37) 

D 



With Lemmas 4 and 5, the NACF and the correlation coefficient can be calculated according to 



h{i) = 






and 



Piii) = 



nil) - {Wi\Y 

E[Xf]-{E[X,]}2' 



(38) 



(39) 



by inserting (36) and (35) into (38), and (36), (35) and (37) into (39), respectively. 

In general, it seems difficult to further simplify (36), (35) and (37). However, we note that 



ln(l + ijjx) 



(40) 



X, w — > 0, 

ln(a'a;), uj — > oo. 

Using (40), we obtain asymptotic closed-form expressions for the NACF, fi(i), and the correlation coefficient, 
pi(i), at low- and high-SNR regimes, as follows. 

1) The Low-SNR Regime: If ry -^ 0, based on (40), (34) can be approximated by 



M 



'^l~2^ -TT'^r.Al), 



rn— 1 



(41) 



which is the same as the low-SNR approximation of X; in a MIMO system with orthogonal space-time block code 



M 



(OSTBC) transmission [20], due to Em=i^™(0 = tr 



h,h; 



= 1 1 H/ 1 1 p. Therefore, the NACF and correlation 



coefficient of interest are equal to those derived for the OSTBC-MIMO system at low SNRs, as stated in the 
following proposition. 

Proposition 1: At the low-SNR regime, the NACF and the correlation coefficient are given by [20] 



h{i) 



NrNt 



1 ' 



NrNt- 

Pi{i) ~ oi- 
2) The High-SNR Regime: If 77 — > 00, based on (40), (34) can be approximated by 

V 



M 



El- 



Nt 



\m{l) 



(42) 
(43) 

(44) 



whose NACF and correlation coefficient are presented in the following theorem. 

Theorem 3: At high SNRs, the NACF and the correlation coefficient are given by (45) and (46), respectively 



fi{i) 



Y^M-l J\f!(m+i/)!ef ^ "^^iF3(M-m.M-rn.M+l,l;AI-7n+l,AI-m+l,N+l;g^) 



- (Em=0 V'J 



=n tPN-m 



^1-1^ 



Ei=oC(2,A^-m) 






^^Inii 



Pii^) 



E™=o '''fM-l^!£l7' 4-^3(M - m, M - m, A/ + 1, 1; M - m + 1, M - m + 1, TV + 1; gf ) 



E™=oC(2,iV-m) 



(45) 



(46) 



where pFg (ai, • • • , Up-, 61, • • • , bq-, z) is the generalized hypergeometric function [24, pp. 1071, 9.14.1], (■(•, •) is the 
Riemann zeta function, given by C,{z, q) — Efc^o ( A\^ [24, pp. 1101, 9.521.1], and ^k is the digamma function 
[24, pp. 954, 8.365.4]. 
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Proof: See Appendix V. D 

Theorem 3 includes the high-SNR approximation for the OSTBC-MIMO system in [20] as a special case. In 
fact, with M = 1, (45) and (46) simplify to the corresponding resutls in [20] by replacing N with MN, i.e., 



j^3i^2(l, 1, 1; 2, MN+\- g^ + kfJV + ln ^) 
^W« z ^—2 -, (47) 

^(^) C(M^AO ' ^^^^ 

where the identity 4F3 (1, 1, 1, 2; 2, 2, MA^ + 1; x) = 3F2 (1, 1, 1; 2, MiV + 1; a;) is used. 

Based on Theorem 3, we conclude that if i^ = and M —> 00, (46) reduces to 

- In (1 - gf) 

lim Pi(*) = T^ ^ ^p \ - S,^o, (49) 

A/-.00 hmp_,oo Lfe=i fe 

where we the first "=" is obtained by collecting the terms in (46), and the second "=" is due to gi < l,i y^O. We 

conjecture that the second "=" of (49) holds for any finite ly at high SNRs, i.e., linij\/_,oo Pj(i) = '^i.o. Vi^ < 00. 

It implies that MIMO IMI is asymptotically uncorrected at high SNRs, if the difference between the numbers of 

Tx and Rx antennas is finite. 

To better understand Theorem 3, the Taylor expansion of (46) and the maximum difference between (43) and 

(46) is listed in Table I, for different values of M and TV. From Table I, the following observations can be made. 

« If 1^ = N — M is fixed, the maximum difference between the low- and high-SNR approximations increases 

when M increases, which is supported by the first four rows of Table I, i.e., (M, N) ~ (1, 1), (2, 2), (3, 3), 

and (4,4). 

> From the last several rows of Table I, i.e., {M,N) = (4,4), (4,8), (4,12) and (4,16), one may conclude 

that if M is fixed, the maximum difference between the low- and high-SNR approximations decreases as ly 

increases. Furthermore, pi{i) can be well approximated by pf , with negligible error for any SNR, when jj is 

not small. 

B. LCR andAOD of MIMO IMI 

The technique developed in Subsection III-B is also valid for calculating the LCR and AOD of MIMO IMI, i.e., 
we can obtain them by replacing (f>\{\th) and ip\{\th.) with (kihh) = Prill > hh} and (^(/th) = Pr{Ii > 
It\i,Ii-i > Ith} in (25) and (31), respectively. Therefore, we only need to calculate, <f>i(Ith) and (^(/th), which 
are presented in the following theorem. 

Theorem 4: At any given level /th, <fe(/th) and (/7r(^th) can be expressed in terms of multiple integrals, given 
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by (50) and (51), respectively. 

M 






^ '^ni j m—1 

p(xi,X2,--- ,Xm) 



2 A/ 



<^Uth) = / " ' / Tv^^^^^ 1 \ 11 ^^^mrfyn 



n:., (i+^)>e.p(..) ^'^! n:::jo^ -k- + .)!.r-^^(i - .?)- cxp (^-^ 



p(a;i,a;2,--- jXa/ ,ai ,i/2,--- ,yM) 

(51) 
froo/- Let {X^}*^^i and {Y^}^,=i be Af unordered eigenvalues of U{l)U'i{l) and H(Z - l)Ht(/ - 1), 

respectively. Then the joint PDF of {Xm}^f=i is given by p(xi,X2, • ■ ■ ,xm) in (50) [30], and the joint PDF of 
{Xjn}fn=i and {Y,n}fn^i is given by p{xi,X2, ■ ■ ■ , XM,yi,y2, ■ • ■ , Vm) in (51) [23]. Moreover, according to (34), 
the event {X; > /th} is equivalent to < nm=i ( ^ + ''jy '" ) > e^"* k which leads to (50). Similarly, it is straightfor- 
ward to see that the two events {Ii > Ith,Ii-i > /th} and {n™=i (l + ^) > e^'N 111=1 (l + ^) > e'"' }, 
have the same probabiUty, which results in (51). D 

Although (50) and (51) can be used to calculate the LCR and AOD of MIMO IMI for small M's, e.g., M = 2, 
via numerical multiple integrals, it is impractical for large A/'s. Fortunately, we can approximate X; as a Gaussian 
random variable for large Af 's and A^'s, which is summarized in the following proposition. 

Proposition 2: If M and N are large, X; can be approximated as a Gaussian random variable with mean /ij = 
E [Ii] and variance o"! = E [if] - {E [X,]}^, where E [X;] and E [if] are given by (35) and (36), respectively [16]- 
[18]. Moreover, we approximate X/ and X;_i by a bivariate Gaussian random vector with mean (E [X;] ,E[X;]) 
and the covariance matrix Ijz ^ aj ( ,.. '^-[ J, where /9r(i) is presented in (39). 

Based on Proposition 2, we have the following theorem for the LCR and AOD of MIMO IMI. 

Theorem 5: Using the Gaussian approximation, we can express the LCR and AOD of MIMO IMI as 

Ai(/th) = 4t /' . ., „. exp(--!%-)d0, (52) 



1 n 

^217, + "- 


/ 

2 \ 


2 sin^ 6 


i-0(/th) 







where /th = "'^~^ is the normalized threshold, and Q(x) — — i= /^ e~~dt is the Gaussian Q-function. 

Proof: See Appendix VI. D 

Theorem 5 requires /ix, ctj and px{l), which can be obtained from (35), (36) and (39). However, for low and 

high SNRs, we may use their corresponding approximations. For high SNRs, they are given by (79), (80) and (46), 

whereas for low SNRs we have /ij = V^R' ^i = ''jy " [20], and f>i{l) ~ gf, obtained from (43). In practice, the 

LCR and AOD at /ij, the ergodic capacity, are of interest, which simplify Theorem 5 considerably. 



Corollary 1: The LCR and AOD of MIMO IMI at the level fix are, respectively, given by 

NM ^ ---^"[/^W] , (54) 

AttTs 

Um) = , ^""^^ ,,,y (55) 

n — 2arcsin[/)r(l)J 

V. Numerical Results and Discussion 

In this paper, a generic power spectrum [19, (8)] [20] is used to simulate time-varying Rayleigh flat fading 
channels with nonisotropic scattering, according to the spectral method [31]. Similar to [20], to verify the accuracy 
of the derived formulas, we consider two types of scattering environments: isotropic scattering and nonisotropic 
scattering with three clusters of scatterers. For nonisotropic scattering, parameters of the three clusters are given by 
[Pi,Ki,6li] = [i,6,0], [P2,K2,6'2] = [i6,f], and [P^^K^i.e^] = [1,8,^], respectively. In addition, in ah the 
simulations, the maximum Doppler frequency Jd is set to lOHz, and 3] = 757^^ seconds. The AoA distributions 
and the corresponding channel correlation coefficients for the above two scattering environments are plotted in Fig. 
2. 

In the following subsections, Monte Carlo simulations are performed to verify NACF, the correlation coefficient, 
LCR and AFD of e/gen-channels and the MIMO IMI of two MIMO systems in the above two propagation 
environments: one is 4 x 4 and the other is 12 x 3. The NACF and the correlation coefficient bear almost the 
same information. The same comment applies to LCR and AFD. Therefore, we only report the simulation results 
for the correlation coefficient and the LCR, to save space. 

A. Eigen- Channels 

In this subsection, the correlation coefficient and the LCR of eigen-channels are considered for both isotropic 
and nonisotropic scattering environments. 

1) Isotropic Scattering: This is Clarke's model [10], with uniform AoA. The comparison between the simulation 
and theoretical results is given in Fig. 3. 

2) Nonisotropic Scattering: This is a general case, with an arbitrary AoA distribution [19][20]. The comparison 
results are presented in Fig. 4. 

In Figs. 3 and 4, the upper left and right subfigures show the correlation coefficient and the LCR of eigen-channels 
in a 4 X 4 MIMO system, respectively, whereas the lower left and right subfigures show the results in a 12 x 3 
MIMO system. In all figures, "Simu." means simulation. In the correlation coefficient plots, "Theo." means they 
are calculated according to (18), and "(fc = I)" denotes the autocorrelation coefficient, whereas "(fc 7^ I)" indicates 
the cross-correlation coefficient. In the LCR plots, "Theo" indicates that the curve is computed using (25)-(28). 

Based on the plots in Figs. 3 and 4, we can see that the derived analytical formulas perfectly match Monte Carlo 
simulations. 



B. MIMO IMI 

In this subsection, the correlation coefficient and the LCR of MIMO IMI are presented for both isotropic and 
nonisotropic scattering environments at low- and high-SNR regimes. In the simulations and theoretical calculations, 
we set r/ = -20 dB for low SNR, and 77 = 30 dB for high SNR. 

1) Isotropic Scattering: For this case, the comparison results are shown in Fig. 5. 

2) Nonisotropic Scattering: The comparison results regarding nonisotropic scattering are given in Fig. 6. 

In Figs. 5 and 6, the upper three subfigures present the correlation coefficient and the LCR of the MIMO IMI in a 
4x4 system. Specifically, the upper left subfigure shows the correlation coefficient at low- and high-SNR regimes, 
the upper middle subfigure gives the LCR of the MIMO IMI at the low-SNR regime, whereas the upper right gives 
the LCR at the high-SNR regime. In addition, the lower three subfigures present the corresponding results in the 
12x3 system. In the correlation coefficient plots, "Theo. (Low SNR)" corresponds to (43), whereas "Theo. (High 
SNR)" corresponds to (46). In the LCR plots, "Theo." means the values are computed from (52), where we used 
low- and high-SNR approximations for the mean /ix and variance (jj, listed immediately after Theorem 5. 
From Figs. 5 and 6, the following observations can be made. 

• Correlation coefficient: If v = max(A^,7VR) — min(7Vr, A'r) is large compared to M = min(A^, A'r), we 
can approximate the correlation coefficient of MIMO IMI by the squared amplitude of the channel correlation 
coefficient for all SNRs, since low- and hign-SNR approximations are very close to each other (see the results 
for the 12 X 3 system). However, if ly is small compared to M, the gap between the low- and high-SNR 
approximations is large (see the results for the 4x4 system). Therefore, we need to resort to the exact 
formulas in (36), (35), (37) and (39) to calculate the accurate values of the correlation coefficient, for not so 
small or large SNRs. For example, at 7/ = 15 dB, the simulation and exact theoretical curves, as well as low- 
and high-SNR approximations are shown in Fig. 7, for the correlation coefficient of the MIMO IMI in a 4 x 4 
system. 
« LCR: The Gaussian approximation works well at both low and high SNRs in large MIMO systems, e.g., the 
considered 12 x 3 channel. But it is not the case in small MIMO systems, say 4x4, where the Gaussian 
approximation has an obvious deviation from the simulation result at high SNR. This is because the central 
limit theorem does not hold for IMI in small MIMO systems^. For this case, we can numerically compute the 
multiple integrals given in (50) and (51), to calculate the LCR. 
« LCR: Compared Fig. 3 and Fig. 4, we find the LCR of an eigen-channel is not sensitive to the scattering 
environment, which is not the case for the LCR of MIMO IMI. Furthermore, based on Figs. 5 and 6, we 
can see that the IMI in a nonisotropic scattering environment has less fluctuations than that in the isotropic 
scattering scenario. 

^In fact, there are obvious differences between the true PDF and the Gaussian approximation in Fig. 1 of [17] at 77 = 15 dB. Larger deviations 
are also observed at higher SNRs, say, 77 = 20 dB, in Fig. 1 of [17]. 
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VI. Conclusion 

In this paper, closed-form expressions for several key second-order statistics such as the autocorrelation func- 
tion, the correlation coefficient, level crossing rate and average fade/outage duration of e/ge«-channels and the 
instantaneous mutual information (IMI) are derived in MIMO time-varying Rayleigh flat fading channels. 

Simulation and analytical results show that the eigen-modes tend to be spatio-temporally uncorrelated in large 
MIMO systems, and the correlation coefficient of the IMI can be well approximated by the squared amplitude of 
the correlation coefficient of the channel, if the difference between the number of Tx and Rx antennas is much 
larger than the minimum number of Tx and Rx antennas. In addition, we have also observed that the LCR of an 
eigen-mode is less sensitive to the scattering environment than the IMI. 

The analytical expressions, supported by Monte Carlo simulations, provide quantitative information regarding the 
dynamic behavior of MIMO channels. They also serve as useful tools for MIMO system designs. For example, 
one may improve the performance of the feedbacked-IMI-based rate scheduler in a multiuser MIMO system by 
exploiting the temporal correlation of the IMI of each user. 

Appendix I 
Proof of Lemma 1 

Although the mean and second moment of X„i{l) were respectively given by (57) and (58) in [16] via a smart 
indirect method, we calculate them directly using its marginal PDF in (7), as follows. 
Using 8.902.2 [24, pp. 1043], we can rewrite (7) as 

where ' mean the derivative with respect to x. With [24, pp. 1062, 8.971.2] 



^^^^ = (i^^^^nL- {[^M-i(^)] 'lU^') - [LU^)]'LL-iix)} , (56) 



and [24, pp. 1062, 8.971.5] 



(56) further reduces to 



[L-;,{x)]' ^ -L:t\ix) (57) 



Lf,(x)=L:;+^(x)-4lll(x), (58) 



P^""^ = (M + z.-^l^!e- {[^^'^-1 (^)]' " ^^'^' ^^)^*^-^ ^4' ^^^'^ 



where the convention Lj^{x) = 0, m < should be used when it is applicable. 
Using (59), we obtain E[A,„(/)] as 

/>oo 

E[A™(0] = / xp{x)dx, 



M + v, 



where the orthogonaHty of Laguerre polynomials [32, pp. 267, 7.414.3] is used, i.e. 

' (k-i-iA^ 

e-^x''Ll{x)mx) = ^-^Skj. 

The last line results in (13), considering N = AI + v. 

By substituting (58) with k = v + \ into (59) and using (61), we can easily obtain (14). 



(61) 



Appendix II 
Proof of Lemma 2 

A. The case of i = 

For i = 0, the value of E[A^(0] is given in (14). 



B. The case of i ^ 

For i ^ 0, we need the following two lemmas. 

Lemma 6: While j, k and ly are non-negative integers, the value of the integral, 

Ii{j,k,i') = J^ x"'^^ e^^ Lj {x)L'j^{x)dx, is given by 



h{j,k,v)=< 



(2fc+iy+l)(fc+iy)! 
fe! ' 

[min(j,fc)+t/+l]! 
[min(j,fc)]! ' 



\J - k\ = 0, 

|j-fc| = l, 

\j -k\>2. 



(62) 



Proof: Using (58), we have 

L^ix)LUx) = L)+\x)Ll+\x) + L)+l{x)Llt\{x) L^+\x)Llt\{x) - ^+l{x)Ll+\x). (63) 

Substitution of (63) into /i(j, fc, v) results in (62), with the aid of (61) and the convention L^{x) = 0, ?ti < 0. D 
Lemma 7: The joint PDF in (6) can be written in the following equivalent form 

, oo A/-1 



j\k\ 



(j + z/)!(fc + i^) 



ef'-''^'^,mxWMLl[x)Li{y) 



y^i 



oX+y J 



j=M k=0 

where p() is the marginal PDF given by (7). 

Proof: By applying the Hille-Hardy formula [32, pp. 185, (46)] 

x + y 



'P{x)p{y), 



(64) 



±7^LUx)LUy)=^'''^ 



k=0 



{k + i^y. 



\-z 



■ cxp 



1 



/.(^r^l,l^l<l, 



to (6), we can obtain (64) after some algebraic manipulations. 
Using Lemmas 6 and 7, it is straightforward to obtain 

1 M\{M-\)\q\ 



.w 



Af2 (M^v)\(M+v-\)\ 
which reduces to (15), based on (62) and N = M ^ v . 



{h{M,M~\,v)X+N'^ 



(65) 

D 

(66) 
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Appendix III 
Proof of Lemma 3 

A. The case of i ~ 

For i = 0, we need the following proposition. 

Proposition 3: If (xi,X2) are a pair of eigenvalues, randomly selected from {Xm{l)}m=i^ then their joint PDF 
is given by [23] 



pixi,X2)'- 



(a;ia;2)^e-(^i+^^) ^^^ p\q\ 



1^^, 



{[L;{x^)L•^^ix2f~L;ix^)L•^^{x,)L;{x2)L•^,ix,)}. (67) 



Note that (67) is different from (6). By reordering the items, we can rewrite (67) as 



3(xi,a;2) = 



M 



M-1 



p{xi)p{x2)^ 



(xia;2)'^e-(-i+"=) "^""^ jlkl 



M{M 



T^ E E (,+4!(fc + .)! ^^^^"^^^^^"^^^^^^"^^^^^"^^- ^''^ 



j=0 A:=0 



Using Lemma 6 and (68), it is easy to obtain 



M 



-N' 



M-1 M(M-l)' 



(69) 



where S = Ejio ^ Efcio^ ^(j+l')!(I^+^^)!" ■ According to (62), we have 

J\/-l A/-2 

5= ^ (2fc+!y+l)2 + 2 J2 {k + l)ik + j^+l), 



k=0 



fc=0 



(70) 



= MN{M + N -1), 



where the last line is derived based on J2l^g k = ^^^^ [24, pp. 2, 0.12L1] and ^^J^q fc^ := »("+iK2»+i) ^24^ 
pp. 2, 0.12L2]. 

Substitution of (70) into (69) proves the first part of Lemma 3, i.e., i = 0. Note that the same result was derived 
in Lemma A of [16] via an indirect method. 

B. The case of i ^ 

Note that {Xm{l)}rn=i and {Xm{l - i)}m=i ^s unordered eigenvalues of Il{l)Ii^l) and H(/ - i)W{l - i), 
respectively, for i ^ 0. So the bivariate PDF of {\m{l), A„(/ — i)}, m ^ n, is the same as that of {Am(/), Am (' — *)}, 
the latter given in (6). Therefore, rm.,n{i) = i'm,mij),i 7^ 0,Vm, n, where rm^m{i) = N^ + '^1^ 7^ 0' is proved 
in Appendix II. 



According to (34), we have 



Appendix IV 
Proof of Lemma 4 



In 1 



E[I;] = ME 

= M \ hi ( 1 



Nr 



A™(/) 



Nt 



(71) 



x 1 p{x)dx, 



where p{x) is given in (7). Substitution of (7) and (29) into (71) results in (35), with the aid of the following 

integral identity [20, (67)] 

0,1 



/ x'^e-^ \n{l+ujx)dx = GI'U- 



0,0,fc+l, 



(72) 



where G is Meijer's G function [24, pp. 1096, 9.301]. 
Similarly, we have 



E[lf] = MI 



lnMl + ^A„(0 



+ Af(M-l)E 



ln(l + ^A,„(OJlnfl + |-A„(0 



m^n 



= M I In^ M + -^x\ p{x)dx + M{M-l) f f ln(l + ^xA\n(l + ^X2\pixi,X2)dxidx2, 

(73) 
where p{x) and p{xi, X2) are given in (7) and (68), respectively. Substitution of (7), (29) and (68) into (73) leads 
us to (36), upon using (72) and the following integral equality [20, (69)] 

k 



x^e-- ln2(l + c.x)dx = 2e^Lo-^^+'^ ^ f ^ j {-if-^Gi;^ I ^ 



-3, -J, -J 
0,-j-l,-j-l,-?-ly 



(74) 



Appendix V 
Proof of Theorem 3 

First we derive the expressions for the first and second moments of X/ in (44), based on the following lemma. 

Lemma 8: Let X = (a;„i_„) be a random matrix with M rows and A^ columns, M < N, where each element is 
a zero mean unit variance complex Gaussian random variable and all the N columns are i.i.d A/-variate random 
vectors with the same M x M positive definite co variance matrix S. The mean and variance of Indct (XX^j are 



E[lndct(XXl' 



M-l 

E 

m=0 



^. 



'N-7n 



IndetE, 



TO 



(75) 
(76) 



Af-l 

var [Indet (XX^)] = ^ ({2, N 

Proof: According to Theorem 1.1 of [33], 2-J'Jdctv has the same distribution as the product of M independent 
X^ random variables with 2A^, 2(A^ —1), • • • , 2(A^ — M + 1) degrees of freedom, respectively. Therefore, we can 



express Indet (XX^) as 



M-l 



Indet (XXt) = X! 1^ ( '^^^^ietS 



(77) 



where the notation = indicates "equal to in distribution", 2y„i is a x^ random variable with 2(A^ — m) degrees of 



, A/-1 



freedom, and {ym},„=o ^^^ independent. 

Based on the results in [20], we have E [In j/„i] = _HAr_,„_i — C and var [In j/,„] — ({2, N — m), where Hk is the 
fc* harmonic number [34, pp. 29, (2.13)], defined by Hk = Y^=i } for fc > 1 with Hq = 0, and C = 0.577215 • • • 



is the Euler-Mascheroni constant [24, pp. xxx]. This completes the proof if we note ipk+i = Hk — C [24, pp. 952, 
8.365.4]. 

It is interesting to observe that the correlation matrix E affects the mean of Indet (XX^) in (75), but has no 
impact on its variance in (76). D 

According to Theorem 1.1.2 of [35] we have 

^ ^^ \det (UiUl),Nn<NT, 

n^™W = i ) ( (78) 

m=i [det (h|H/ j , Nr > Nt. 

By applying Lemma 8 and (78) to (44) with S = Im, it is straightforward to write the mean and variance of X; as 

E[Ii]^J2 "I'N-rn + M\n^, (79) 

and 

M-l 

var[Ii]« Y, C%N-m), (80) 

?n— 

respectively. These two are consistent with the results in [16], where an implicit complex extension of Theorem 
3.3.4 of [35] was used. Clearly, the second moment of Ii is given by 

M-l /M-l \^ 

E[lf]« ^ C(2,iV-™)+^7^jv-™ + A/ln^ . (81) 

?n— \rn— / 

For calculating the autocorrelation of X;, we need the following lemma. 

Lemma 9: With j, k and i/ as non-negative integers and j ^ k, the value of the integral 

hij, k, v) = /g°°(lna;) x" e~^ L'^ {x)L'j[{x)dx is given by 

w ■ . N [mm{j,k) + i']l 

^-(^'^''-)- [min(,,fe)]![minO-fc)^maxafc)] - ^''^ 

Proof: First we consider j > k. Substitution of L'^{x) with (29) into /2(j, k, v) gives 

Uj^k,.) = Y: (l^_''Y-^r{lnx)xP+^e-'^L^{x)dx, 

where the second line comes from 2.19.6.2 [36, pp. 469]. Using X;^o ("«') = CVli^) [^4, pp. 4, 0.151.1], we 
have E^=o ^^'^='^10 i'tl-'i')^^ - k - l)\ = (/J(j - k 1)! = ,^, which reduces (83) to 



Similarly, for j < k, we obtain 



h{j.k,v)= .;^ ; (85) 

jKj - k) 



Combination of (84) and (85) results in (82). D 
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Now we proceed to prove (45) and (46). Based on the high-SNR approximation of X; in (44), we have 

M M 



m — 1 n— 1 



' ■qX^.Q.) , r]X„{l-i) 



Nt 



Nt 



^^ ??Ai(0 ^^^7]Xi{l-i) 



Nt Nt 



(86) 



^ Rr 



ln2 ^ + 21n-J^E[lnA(0] +nnA(*) 

I\t Jvt 



where rinA(i) = E [In A(/) In A(Z — i)]. Using (64) and Lemma 9, r\nx{i) can be evaluated as 
By substituting (82) and (87) into (86), we obtain 

M-l 



(87) 



(88) 



fc=0 



where S{k, v, g.) = E,^,/ ilgtS^fc)^ ' ^"^^ ^t^'l ^ ^^^ 



In 






is approximated by (79). By introducing a 



new variable p = j — AI in S{k, v^ qi) and using the Pochhammer symbol (a;)„ = x{x + 1) • • • (x + 71 — 1), we can 
rewrite S{k,u, Qi) as 

.2p 



S{k,v,Qi) 



M\{k + J.)!ef '"-'=) ^ [(M - fc)p]2 (M + l)p(l)p e, 



■E' 



(M - kYN\k\ ^ [(M _ fc + 1)J' [N + l)p p! 
M!(fc + :/)!gf *'-'=) 



(89) 



-4i^3(Af-fc,Af-fc,A/ + l,l;A/-fc+l,Af-fc + l,iV+l;£i-), 



(A/ - /c)2Ar!fc! 

where N — Af + v, and the last line comes from the definition of the generalized hypergeometric function [24, pp. 
1071, 9.14.1]. 

Substitution of (81), (88) and (89) into (38) results in (45). Similarly, with (79), (80), (88) and (89), (39) reduces 
to (46). 



Appendix VI 
Proof of Theorem 5 

To simplify the notation, we set X = X;, Y = X;_i, and p = pi{l). According to Proposition 2, we have the 
PDF of X and the joint PDF of X and Y as 



p{x) 



I-KOI 



exp 



{x - Mr) 
2a| 



and 



Pix,y) 



exp 


(a:-Aix)^ + (j/-Mx)^-2p(2;-A'x)(!/-A'z) 


[ 2a|(l-p2) J 



(90) 



(91) 



V27r(l-p2)oi 

In what follows, we calculate (j)i{Ith) ~ Jj p{x)dx and ipi{Ith) = Jj Jj p{x,y)dxdy for the cases of /th > Mr 
and /th < Mr- 
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A. The Case of /th > A*r 

According to (4.2) [26] we obtain 



^(/th) '"•-"" i /^ cxpf --^ I d0, (92) 

TT Jo \ 2 sm 61 / 



where /th — "* ^^ . Similarly, using (4.18) [26] and the following equality 



oi 



I- pj 4 
we obtain 



, ,l + p\ TT arcsin (p) ,„^^ 

arctan J— ^ =- + — ^^, (93) 



C3in(p) 



^(/,,)'"^^^^i/"^ ^ exp(--i%-)d0. (94) 

TT Jo \ 2 sni 9j 

Substitution of (92) and (94) into (25) results in (52). Moreover, ii(/tii) = /_^'^ p{x)dx = 1 - Q f/th) ■ By plugging 
Fi{Itb) and (52) into (31) we obtain (53). 

B. The Case of /th < Mr 

For this case, using the results in (92) and (94) and the symmetry of the Gaussian PDF, i.e., the integral equality 
f —^e~~dt ^ f -^=e~~dt, it is straightforward to obtain 

<fe(/th) '"='^" 1 - - /^ °^pf-^7^1 ^^' (95) 

and 



<^(/th)'-i^^l--/ exp--i%- d0+i/ exp--i%-d0. (96) 

We obtain (52) by substituting (95) and (96) into (25). Similarly, we get (53) easily by plugging Fx{Ith) and (52) 
into (31). 
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TABLE I 

Taylor Expansion of (46) and the Maximum Difference between 

(43) AND (46) FOR Different (M, A'')'s 



(M, N) 


Taylor Series of (46) 


maxo<e.<i|(43)-(46)|t 


(1,1) 


0.608^2 +0.152e| + O(ef) 


0.160 


(2,2) 


0.437^2+0.218^1 +0(ef) 


0.230 


(3,3) 


0.372^2+0.186^1 +0(ef) 


0.274 


(4,4) 


0.337^2+0.168^4 + 0(ef) 


0.304 


(4,8) 


0.725^2 +0.178^4 + 0(^6) 


0.085 








(4, 12) 


0.824^2+0.135^1 + 0(^6) 


0.050 








(4, 16) 


0.870^2+0.1074 + 0(^6) 


0.036 



■f The maximum difference is calculated via the function FindMaximum 
in Mathematica® . 
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Fig. 1. (a) A MIMO channel with Nt transmit and Nr receive antennas; (6) The equivalent M parallel SISO channel representation. 
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Fig. 2. The AoA distributions for two scattering examples and the corresponding channel correlation coefficients. 
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Fig. 3. The correlation coefficient and the LCR of an e(ge«-channel, in 4 X 4 and 12 X 3 MIMO systems with isotropic scattering. 
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4x4 Nonisotropic 
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Fig. 4. The correlation coefficient and the LCR of an ei'gen-channel, in 4 X 4 and 12 X 3 MIMO systems with nonisotropic scattering. 
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Fig. 5. The correlation coefficient and the LCR of the MIMO IMI at low- and high-SNR regimes, in 4 X 4 and 12 X 3 MIMO systems with 
isotropic scattering. 
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Fig. 6. The correlation coefficient and the LCR of the MIMO IMI at low- and high-SNR regimes, in 4 X 4 and 12 X 3 MIMO systems with 
nonisotropic scattering. 
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Fig. 7. The correlation coefficient of the MIMO IMI at 77 = 15 dB (moderate SNR), in a 4 X 4 system with isotropic scattering. 



